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Aims. The aim is to elucidate the astrophysical conditions required for generating large amounts of dust in massive starburst galaxies 
at high redshift. 

Methods. We have developed a numerical galactic chemical evolution model. The model is constructed such that the effect of a wide 
range of parameters can be investigated. It takes into account results from stellar evolution models, a differentiation between diverse 
types of core collapse supernovae (CCSN), and the contribution of asymptotic giant branch (AGB) stars in the mass range 3-8 M Q . 
We consider the lifetime-dependent yield injection into the interstellar medium (ISM) by all sources, and dust destruction due to 
supernova (SN) shocks in the ISM. We ascertain the temporal progression of the dust mass and the dust-to-gas and dust-to-metal mass 
ratios, as well as other physical properties of a galaxy, and study their dependence on the mass of the galaxy, the initial mass function 
(IMF), dust production efficiencies, and dust destruction in the ISM. 

Results. The amount of dust and the physical properties of a galaxy strongly depend on the initial gas mass available. Overall, while 
the total amount of dust produced increases with galaxy mass, the detailed outcome depends on the SN dust production efficiency, 
the IMF, and the strength of dust destruction in the ISM. Dust masses are higher for IMFs biased towards higher stellar masses, even 
though these IMFs are more strongly affected by dust destruction in the ISM. The sensitivity to the IMF increases as the mass of the 
galaxy decreases. SNe are primarily responsible for a significant enrichment with dust at early epochs (< 200 Myr). Dust production 
with a dominant contribution by AGB stars is found to be insufficient to account for dust masses in excess of 10 s M G within 400 Myr 
after starburst. 

Conclusions. We find that galaxies with initial gas masses between 1-5 x 10" M Q are massive enough to enable production of dust 
masses > 10 8 M Q . Our preferred scenario is dominated by SN dust production in combination with top-heavy IMFs and moderate 
dust destruction in the ISM. 
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1. Introduction 

Modelling the evolution of dust in galaxies is a key ingredient 
in understanding the origin of the high observed dust masses in 
high-redshift galaxies and quasars (QSOs). Dust masses > 10 8 
M have been derived from observations in QSOs at redshift z 

> 6 (e g.JBertoldi et al]l2003at iRobson et al|2004l:|Beelen et alJ 
2006; iHines et alj|200a iMichalowski et aljboiObl) along with 
high star formation rates (SFR) up to a few times 10 2 ~ 3 M yr 1 
(e.g. JWalter et al]l2004r. IWang 'et al.ll2010h . Additionally, QSOs 
at z > 6 harbou r supermassive black holes (SMBHs) with masses 

> 10 9 I vL (e.g-JWillott et alj2003UVestergaardl2004HJiang et all 
2006). Kawak atu & Wadal (120081 l2009h showed that to form a 
SMBH > 10 9 M Q , a high mass supply of > 10 10 " 1 1 M Q is needed. 
These requirements, together with derived molecular gas masses 
from CO line emission measurements i n excess of 10 10 M R 
(e.g.. ICox et al.ll2002t ICarilli et al ll200l iBertoldi et alJl2003bt 
Walt er et al J 12003b iRiechers et alJl200sTlWang et alJl2010h . set 
significant constraints on the physical properties of the host 
galaxies. These in turn have implications for the origin and evo- 
lution of dust. Furthermore, a tendency toward increased dust 
attenuation with hi gher galaxy masses for syst ems at z ~ 6-8 
has been found bv lSchaerer & de Barrosl d201Q|) . 

Analytical a nd numerica l mode ls for dust evolution have 
been developed. iDwek i t al. (2007) propose 1 M Q of dust per 



SN will be necessary to account for dust masses in high-z QSOs, 
while contemplating SNe as the only source. Such high dust 
masses for SNe contradict derived dust masses from nearby 
SNe and SN remnants, w hich on average reveal a few times 
10-4 10-2 Mo f dust (e.g . JWooden et alJl993t lElmhamdi et all 
[20031: iT^mim et alj|200a iMeikle et all 120071: iRho et alJ 12 008: 
kotak et alj|2009t ISibthorpe et alj|2010t iBarlow et alj|2010h . A 
review of observationally and theoretically derived dust from 
stellar sources is provided by Gall et al. (in prep, herafter 
GAH11). 

The issue of whether SNe produce large amounts of dust 
is unclear. Other sou rces of dust such as AGB stars have been 
taken into account (Morgan & Edmunds 2003; Valiante ~et al.l 
| 2009h in chemical e volution models of high-redshift galaxies. 
IValiante etall d2009l) claim that with the contribution of AGB 
stars, 10 8 M Q of dust can be reached, with AGB stars dom- 
inating the dust production. AGB stars (0.85-8 M ) are the 
main source of dust in the present universe, but only stars 
with masses > 3 M Q are likely to contribute at z > 6 (e.g., 
Marchenko 2006). Evidence that metal-deficient AGB stars 
also undergo strong mass loss and are abl e to efficiently pro- 
duce dust is supported ob s ervationally (e . g.JZijlstra et al.ll 2006: 
Groenewegen et al. 2007J; lLagadec et al.l 120071 : Matsuura et al, 
2007; ISloan et al.ll2009h and theoretically (e.g.. IWachter et al 
2008i iMattsson et al.l 120081) . However, the theoretical mod- 
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els dDwek et al.ll2007t lMorgan & Edmund sll2003l: IValiante et al l 
2009) greatly differ with respect to the assumptions made for the 
mass of the galaxy and dust contribution from stellar sources, as 
well as the treatment of the star formation. Thus the origin of 
dust and its evolution remain unclear. 

In GAH11 we discuss plausible dust production efficiency 
limits for stellar sources between 3 M and 40 M , and deter- 
mined the dust productivity of these sources for a single stellar 
population. In this paper, we investigate the evolution of dust 
in high-z galaxies. We develop a numerical chemical evolution 
model, which allows exploration of the physical parameter space 
for galaxies at z > 5-6. Different types of core collapse super- 
novae are specified and the contribution from AGB stars and the 
impact of SMBHs are taken into account. We furthermore fol- 
low the evolution of some physical properties of these galaxies. 
The main parameters varied in the model are the IMF, the mass 
of the galaxy, yields for SNe, and the strength of dust destruc- 
tion in the ISM, as well as the dust production efficiency limits. 
Models with or without the SMBH formation are considered. 

The paper is arranged as follows. In Sect. [2] the equations 
used to construct the model are developed. We discuss the model 
parameters and their possible values in Sect. [3] A detailed anal- 
ysis of the results is presented in Sect. |4] which is followed by a 
discussion in Sect.[5]and our conclusions of this work in Sect. [6] 

2. Modelling the evolution of dust in starburst 
galaxies 

In this section we formulate the equations needed to follow 
a galaxy's time-dependent evolution in a self-consiste nt nu- 
meric al model. The main basic logic is adopted from iTinslevi 
(1980, and references therein), which has also been used in 
other chemical evolution models to study du st in galaxies (e.g., 
Morgan & Edmunds 20031 iDwek et aDl2007l) . We focus on an 
elaborate treatment of dust from different types of CCSNe 
and AGB stars. In order to calculate the amount of dust 
from these sources we use the dust production efficiencies de- 
scribed in Sect. 13.41 In particular, the lifetime-dependent de- 
layed dust and gas injection from AGB stars and SNe is taken 
into account. The metallicity-depe n dent lifetimes of all stars 
are taken from Schaller et al ] d!992l) . ISchaerer et alJ (I 19931) . and 
ICharbonnel et all (119931) . In GAH1 1 we show that the variation 
of the lifetime with metallicity is minimal. We therefore calcu- 
late and use a metallicity-averaged lifetime for all stars. The re- 
cycled gaseous material is defined as the remaining ejected stel- 
lar yields from all massive stars in the mass range 3-100 M , 
which has not been incorporated into dust grains. This also in- 
cludes the stellar feedback from very massive stars (> 30^40 
M ) in the form of stellar winds. Dust and gas are assumed to be 
released instantaneously after the death of the stars. 

We strictly treat the elements in the gas and solid phases 
separately, while paying attention to their interplay. Thus, we 
define M g (f) as the total mass of elements in the ISM, which 
are in the gas phase and M&{t) as the total amount of elements 
in the solid dust phase. The mass of the ISM is defined as 
Mism(0 = M„(t) + M&(i), which in the literature is often referred 
to as the 'gas mass'. The total amount of dust in our models 
is solely calculated from the dust contributions from SNe and 
AGB stars. Thus, no further growth in the ISM is contemplated. 
However, dust destruction in the ISM through SN shocks is taken 
into account. The formation of the SMBH is considered as a sim- 
ple sink for dust, gas, and metals. 

We assume a so-called 'closed box' model; i.e., the effect of 
infalling and outflowing gas in the galactic system is neglected. 



Although galaxies may not evolve in such a simple manner, this 
assumption is plausible since massive starburst galaxies are as- 
sumed to have SFRs iff(t) > 10 3 M Q yr _1 . Infall of neutral gas 
would only affect the system when the infall rate is compara- 
ble to the SFR. In this case a large gas reservoir needs to be 
present in the vicinity of the galaxy already. Besides, infall rates 
for high-z galaxies are not known. We further assume the ISM 
to be homogeneously mixed, and we evolve our model only up 
to the first Gyr. Examples of 'closed box' models being suffi- 
ciently accurate in clude the works of iFraver et alJ (1 1999b and 
iTecza et all d2004l) for the luminous and massive submillimetre 
galaxy SMMJ1401 1+0252 at z = 2.565. 

2. 1. Basic considerations 

The initial mass functions (IMF) <p(m) is normalized to unity in 
the mass interval [m\, m^X as 



J in] 



m <p(m) dm — 1 , 



(1) 



where m\ and mi are the lower and upper limits of the IMF (see 
Sect. EH*. 

The relation between the total m ass Mism(0 and the SFR 
tKO is given by the Kennicutt law dKennicuttl 1 1998b . where 
(/r(f) oc Mism(0*- Analogously to IDwek et al.l (|2007), we apply 
the following notation to calculate the SFR 



M ISM (f) 



(2) 



where ip m \ is the initial SFR, M m \ the initial gas mass of the 
galaxy, and Mism(0 the mass of the ISM. The value of k is be- 
tween 1 and 2, so we set k = 1.5. This va l ue has often been 
assumed in other models (e.g.. lDwek1ll998t IDwek et al.l 12007b 
ICaluraetal.ll2008l 



2.2. Equations for AGB stars and supernovae 

The amount of dust from all stellar sources released into the ISM 
per unit time is simply calculated as the amount of dust produced 
by the stellar sources times the source rate. The considered dust 
producing stellar sources are AGB stars and CCSNe. We account 
for a potentially diverse dust contribution by different SNe sub- 
types and distinguish between Type IIP SNe and the remaining 
Type II subtypes. Types lb and Ic SNe, collectively referred to 
as Ib/c, are not considered as dust producing SNe. However they 
inject their stellar yields into the ISM. A specification of these 
sources and their lower and upper stellar mass limits, and 
mu(i), are discussed in Sect. 13.31 In the following equations these 
sources are indicated by the subscript i = AGB, IIP, II, Ib/c. 
The AGB and SN rate Ri(t) calculates as 



Ri(t) 



= J ' 



if/(t — t) <p(m) dm, 



(3) 



where t = t(to) is the lifetime of a star with a zero-age main 
sequence (ZAMS) mass m; i.e., the star was born at time (t - t) 
when it dies at time t . It is evident that stars only contribute when 
the condition t - r > is fulfilled. 

For SNe, the time of releasing the total produced elements is 
assumed to take place right after explosion. The main sequence 
lifetime for AGB stars is defined as the time until the end of 
the early AGB phase, which can take up to several 100 Myr. 
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However, the most efficient mass loss phase itself is less than 
1 Myr at the very end of the AGB phase and is relatively short 
compared to the total lifetime of AGB stars. Thus, we make the 
same approximation as for SNe: All produced elements and dust 
are released instantaneously after the main sequence lifetime. 

For each kind of source the total produced dust per unit time 
is calculated as 

£ df( f ) = (Jz+Jwind Z(f-r)) eiim, Z) W-t) <p(m) dm, (4) 

where Yz = Yz(m,Z) with Z = Z(t - r) is the mass (m) and 
metallicity (Z) dependent amount of ejected heavy elements per 
star. The metallicity Z is defined as the metallicity with which the 
star was born at a time (t - r). The parameter e\(m, Z) is defined 
as the dust production efficiency. The assumed efficiencies are 
described further in Sect. [14] For Types Ib/c SNe, E^Jt) = 0. 

The total amount of mass lost in stellar winds prior to ex- 
plosion, Iwind, caused by mass loss during stellar evolution of 
SNe, is calculated as 

JVind = m - M fin , (5) 

where Jwind = Jwind^Z) and Mf\ n = Mf\ n (m,Z). The latter 
is the final mass of a SN before explosion. Type II SN suffer 
strong mass loss leading to the formation of dense circumstellar 
disks. Dust has been found in these disks for some SNe (GAH1 1 
and references therein). However, the actual amount of metals 
in these disks is not known and data are not available, thus we 
cannot account for them. In the case of AGB stars and for SN 
models where no mass loss prescription is available Fwind = 0. 

Besides the elements bound in dust grains, elements in the 
gas phase are also released into the ISM. The produced amount 
of elements in gaseous form is calculated as 

(Ye + YwiNoMt ~ t) 0(m) dm - E%™(t), (6) 

m m 

where Ye = Ye{iii,Z) is the amount of all ejected elements per 
star. 

The total mass of heavy elements released into the ISM per 
unit time is calculated as 

Jr*n m 
(Yz + Iwind Z(? - t)) ifr(t - t) <f>(m) dm. (7) 
»lL(i) 

We include the possibility of dust destruction in the SN remnant 
(SNR) due to reverse shock interaction. D ust destruction time 
scales up t o 10 4 yr have been predicted bv lBianchi & Schneider! 
d2007l) and lNozawa et all d2007ll2010t) . However, this timescale 
is relatively short in comparison to the lifetime of a star, so we 
make the approximation that dust is destroyed immediately after 
formation. 

We define the parameter £sn as the SN dust destruction fac- 
tor; i.e., the destroyed mass of dust per unit time of all SNe is 

£ d d f(f) = £ d 7 d (0ftN. (8) 
This term applies only to Type IIP and Type II SNe, thus i = IIP, 

n. 

The final SN dust injection rate per unit time is calculated as 

£d, SN (o= 2 <r d «-£ d d r(o, (9) 

i=np,n 

while the final AGB dust injection rate is 

£d,AGB(0 = <a G b(0. (10) 



2.2.1. Recycled gaseous material 

The recycled material from SNe, and AGB stars consists of all 
the mass of the elements not being incorporated into dust grains, 
thus the material is in the gas phase. Very massive stars ending as 
BHs may contribute with their stellar winds to the recycled ma- 
terial. We refer to stars that directly form a BH as the 'remaining 
stars'. Pertaining to the short lifetime of very massive stars and 
the resulting short duration of the wind phase, we assume that all 
the elements lost in the wind phase are released after the death 
of the star. The total mass of the released elements per unit time 
of the remaining stars is 

J^™U(R) 
X e , K if,(t-T)<f>(m)dm, (11) 
«1L(R) 

where X„r = X„^(m,Z) is the mass of elements released into 
the ISM per star, while the subscript 'R' stands for 'remaining 
stars' . 

For these remaining stars the following tw o scenarios are 
possib le: (1) a supernova without display (Eldridge & Tout! 
2004) occurs even if a BH is formed and elements are ejected or 
(2) no SN occurs because no SN shock is launched. For the first 
case, the term X gt R can either be substituted with (Ye + ^windX 
or if Ye and Fwind are not known, the common approximation 
°f ^g,R = m-M Km can be made. The mass M rem = M mm (m, Z) is 
the remnant mass of a star. In the second case we assume that no 
nucleosynthesis products will be ejected. A contribution to the 
gas household in the ISM comes solely from the stellar winds, 
thus X s r = Fwind- 

The mass of released metals from these stars per unit time is 
given by 

Xz,Ri[f(t-T)<f>(m)dm, (12) 

l«L(R) 

where Xz,r = Xz,R(m,Z) is the mass of the ejected heavy ele- 
ments per star. For the first case, Xz,r = Yz + ^wind Z(f - r). For 
the second, Xz,r = ^wind Z(f - t). 

By extending the subscript i to i = AGB, IIP, II, Ib/c, R, the 
total returned mass of gaseous material to the ISM per unit time 
calculates as 

£ g (0 = 2 £ g.i(')+ J] £ d d f(0- (13) 

i i=IIP,U 

It is important to note that elements locked up in dust grains will 
return to the gas phase when dust destruction takes place. As a 
result, the amount of destroyed dust in SNRs must be added to 
the total amount of gas in the ISM (second term in Eq.[T3]>. 

2.3. The evolution of dust and gas in the galaxy 

The chemical evolution of a galaxy is mainly determined by the 
equations balancing the net amount of gas, dust, and heavy ele- 
ments. 

2.3.1 . Effect of a super massive black hole 

We take the effect of a SMBH into account as an additional 
sink for the gas, dust, and heavy elements. We base our as- 
sump tions for the treatment of the BH on the theoretical work 
of iKawakatu & Wadal (l2008l l2009h . Super Eddington growth is 
required to form a SMBH. This necessitates a large mass sup- 
ply of ~ 10 10 ~ n M on a short supply timescale of fsMBHsup ~ 
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10 8 yr from the host gala xy to a massive circumnuclear disk. 
iKawakatu & Wadal (l2009h find that the final SMBH mass is 
around 1-10 % of the supply mass M B Hsup- 

We do not treat the formation of the disk in detail, but assume 
that the supply mass needed is equal to the initial mass M m i of 
our considered systems. We therefore only take the overall mass 
of the ISM needed to form the SMBH into account. A simple 
constant growth rate is calculated as 



SMBH - 



MsMBH 
fSMBH 



where Msmbh is the mass of the final SMBH, and *Psmbh = 
const, for t < f SM BH, whereas ^smbh = for f > f SM BH- The 
growth timescale to build up the SMBH is equal to the supply 
timescale ?smbh = ?smbhsup- The onset of the SMBH formation 
coincides with the onset of starburst of the whole galaxy. Where 
the SMBH is not taken into account, the growth rate is set to zero 

(^SMBH = 0). 



2.3.2. The amount of dust in the ISM 

The evolution of the amount of dust Md(f) in the galaxy is 



dM d (t) 
dt 



£d,SN(0 + £d,AGB(0 - E D (t). 



(15) 



The first and second terms are the dust injection rates from SNe 
and AGB stars contributing to the increase in the dust household 
in the ISM. The third term Eu(t) is defined as the total dust de- 
struction rate. It determines the dust reduction through astration, 
as well as through the SMBH formation and destruction in the 
ISM caused by SN shocks, if considered. The total dust destruc- 
tion rate is calculated as 



E D (t) = T] d (t) (t/Kr) + Tsmbh + M a R 5 n(t)). 



(16) 



The variable Tja(t) = Md(f)/.MiSM(0 is the fraction of dust in the 
ISM and will be referred to as the 'dust-to-gas mass ratio'. For 
simplicity we make the assumption that the produced dust will 
be immediately mixed with the material in the ISM. The dust 
destruction in the ISM through SN shocks will be discussed in 
more detail in Sect. 13.51 

2.3.3. The amount of gas in the ISM 

The temporal evolution of the gas content in the galaxy is calcu- 
lated as 



dM g (f) 
dt 



= £ g (f) + 77d(OM clJ R SN (f) 
-(I-^OMO + ^smbh). 



(17) 



Here, E g (t) is the recycled gaseous material (see Eq. [TJJ. The 
second term is the amount of destroyed dust in gaseous form 
(see Sect. l2.3.2i i. The third term accounts for the depletion of the 
gas in the ISM through incorporation into stars and the loss to the 
SMBH. The fraction of gas in the ISM is expressed as (1 -J]i(t)). 

2.3.4. Metallicity 

Owing to the separation of the ISM mass into the material locked 
in dust and gas, respectively, it is necessary to formally take care 
of the transition of elements from the dust phase to the gas phase 



(due to destruction). For the evolution of heavy elements we do 
not distinguish between their chemical states. 

The equation for the evolution of the total amount of heavy 
elements is formulated as 



dM z (t) 
dt 

where 



^zW-^zOOOKO + ^smbh), 



(14) £z(f) = J] £z - l(f) 



(18) 



(19) 



is the total ejected mass of heavy elements per unit time from 
all considered sources. The last term in Eq. [T8l determines the 
reduction of heavy elements due to astration and the loss to the 
SMBH. 

The total metallicity of the system is defined as Z(f) 
if/it) = M z (t)/M 1SM (t). The fraction of metals in the ISM which 
are bound in dust grains is calculated as rjzd(t) = Md(?)/Mz(?)- 
The amount of metals in the gas phase is 

M z , g (*) = (l-T]rt(t))Mz(t). (20) 

The gas phase metallicity is given as r/z g (t) = M z>g (f)/M g (f). 

3. Model parameters 

In this section we describe the prime model parameters, along 
with the values used in this study. In particular, we consider the 
initial conditions of the galaxy, the IMF, stellar yields, and the 
destruction rates of dust in the ISM. These characterize the sys- 
tem and significantly influence the evolution of gas, dust, and 
metals. In addition we define some switches, which specify var- 
ious possibilities for some model parameters. All parameters 
and their considered values, as well as the possibilities for the 
switches, are listed in Table[2] 

3. 1. Initial conditions 

The model is defined by the initial values of the parameters. It 
neither depends on nor is influenced by additional input from 
other models; i.e., it does not depend on a cosmological model. 
Therefore it can be applied to any galaxy within the accuracy 
limit of a 'closed-box' treatment. We are mainly interested in 
massive high-redshift galaxies in which high dust masses, stellar 
masses, and H2 masses have been inferred from observations. 
The parameters of our computed models are therefore tuned to 
such galaxies. 

One of the main parameters is the baryonic initial gas mass 
Mi n i, which is equal to the total mass of the galaxy in baryons. A 
relation between M; n j and the mass of the dark matter halo Mdm 
hosting such systems is given through M; n j - Q.t,/Q. m Mom- In 
this work we consider four different massive galaxies with M m \ 
= 1.3 x 10 12 M , M ini = 5 x 10 11 M , M ini = 1 x 10 11 M Q , and 
Mni = 5 x 10 10 M . 

In GAH11 we argue that the very first population of stars 
(so-called Pop III stars) are not likely to be the main sources of 
high dust masses at high redshift. Thus we consider only the next 
generations of stars (Pops II or I). The formation of these stars 
takes place as soon as a critical metallicity of Z„ ~ 10~ 6 -10~ 4 Z p 
( Bro mm & Loebll2003t ISchneider et aT]l2006t lTumlinsonll200 



is reached in the star-forming region. In this regard we assume 
an initial metallicity in accordance with the critical metallicity 
ofZi„i=Z cr = 10- 6 Z o . 
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Table 1. IMF and the adopted values 



IMF 


a 


m\ 


m-i 


ra ch 


Salpeter 


2.35 


0.1 


100 




Mass-heavy 


2.35 


1.0 


100 




Top-heavy 


1.5 


0.1 


100 




Larson 1 


1.35 


0.1 


100 


0.35 


Larson 2 


1.35 


0.1 


100 


10.0 



Table 3. Definitions of stellar types 



Type Af H (w, Z) [M g ] M CBK (m, ZjjMgL 



Type II-P > 2 < 1 

Remaining Type II > 0. 1 and < 2 < 1 

Typelb/c <0.1 <1 

Remaining stars — > 1 



Pertaining to the rather high derived star formation rates 
from observations of some high-z mass i ve galaxies and QSOs 
(e.g iFraver et alJ Il999t iBertoldi et ail l2003at iRiechers et all 
120071) . we consider an initial SFR of ^ ini = 1 X 10 3 M Q yr -1 . The 
evolution is determined using the Kennicutt law as described in 
Sect. 12. 11 Eq.[2] We also consider a case of constant SFR where 
iKO = <Aini = 1 x ltfMsyr 1 . 

In our model the onset of starburst is not directly connected 
to redshift, so, the age of the galaxy is identical to the evolu- 
tionary time after starburst. For dusty galaxies seen at redshift 
5-6 the earliest onset of starburst with very high SFRs can be 
considered to have taken place at z — 10. For a ACDM universe 
with H p = 70 km s" 1 Mpc ', Q. A = 0.73 and Cl m = 0.27 and Q. b 
= 0.04 (Spergel et al. 2003), the evolutionary time of interest for 
building up high dust masses possibly lies then within 400-500 
Myr. We have computed all models presented in this paper up to 
an age of the galaxy of f max = 1 Gyr. 

3.2. The initial mass function 

The initial mass function is an important parameter influenc- 
ing the evolution of dust, gas, and metals in a galaxy. It deter- 
mines the mass distribution of a population of stars with a certain 
ZAMS mass. In models of dust evolution in galaxies and h igh-z 
quasars (e.g.. Morg an & Edmunds 2003; D wek et al]l2007l) . an 
IMF first proposed by Salpeter (1955) is often used. However, 
for starburst galaxie s there is also observational evidence of a 



top-heavy IMF (e.g. 



Dabringh ausen et al. 



Doane & Mathewslll993UTumUnsonl l2006 ; 
120091: lHabergh am et al]l2010l) . 



In view of this, we consider a set of five different IMFs. The 
power-law IMFs (Salpeter, mass heavy and top-h eavy) have the 
form <p(m) oc m~ a . The lognormal Larson IMFs dLarson| [l998) 
are given as (p(m) oc m~^ a+1 ' exp(— m C h/m), where m c h is the char- 
acteristic mass. Their parameters are presented in TableQ] 

The influence of the diverse IMFs, particularly on the 
dust production rates from stellar sources, is demonstrated in 
GAH11. 

3.3. Stellar yields 

The model is adapted to published results from stellar evolution 
models. Using these models the metallicity dependent upper and 
lower mass limits of the mass range of diverse SNe subtypes 
can be determined. For comparison to the usually adopted treat- 
ment of SNe we calculate simplified models with a fixed mass 
range for CCSNe between 8^4-0 M and use stellar yields from 
nucleosynthesis calculations for SNe. For SNe we adopt three 
different structural models. 



3.3.1 . SN yields from stellar evolution models 

We adopt the stellar yields for the hydrogen mass Mu(m, Z) in 
the envelope, the He core mass M COK (m, Z), and the final mass 



Mfln(m,Z) of the progenitors, as well as Y^{m,Z) and Yz(m,Z) 
from the single stellar evolution models by Eldridge et al. (2008, 
herafter EIT08I). These yield s are ba sed on previous works 
dEldridge& Toutl I2004 l2005l) IEIT08I adopt a mas s loss pre- 
scriptio n using th e rates of [cteJager et al. (1988) with the 
rates of IVink et all (1200 ll) for pre-Wolf-Rayet (WR) and from 
iNugis & Lame rs (2000) for WR evolution including overshoot- 
ing. In their study, this model is closest to a set of SN progenitor 
obs ervations. 

lEldridge & Toutl (120041 120051) also show that all stars end up 
as II-P SN at metallicities roughly below Z = 10~ 4 , while with 
increasing metallicity the upper mass limit for II-Ps decreases; 
i.e., at solar metallicity it is at roughly 28 M Q . Type II-P SNe 
are per definition all stars that have retained at leas t 2 M of 
hydro gen in their envelopes at their pre SN stage dHeger et al.l 
2003). Since mass loss is more efficient at higher metallicities, 
the upper mass limit for II-P supernovae decreases. 

With increasing metallicity further types of CCSNe such as 
IIL, lb, and Ic SNe arise. The upper limit for IILs is defined by 
the small hydrogen fraction of ~ 0. 1 M Q in the envelope. In case 
little or no hydrogen in the envelope of the progenitor is present, 
SNe appear as Types lb or Ic. The SNe arising from the higher 
mass end of all Type lis are the subtypes lib or Iln. These types 
are difficult to fi t into a quiescent mass loss prescription and are 
not specified by EIT08. Thus we simply assume that these sub- 
types may be part o f the fraction of SNe, wh ich is determined as 
IIL in the models of lEldridge &~fo ut (2004). These SNe will be 
collectively referred to as the remaining Type II supernovae. 

We define the conditions for determining the lower msL,© 
and upper OTsu.m boundarie s of th e different SN subtypes anal- 
ogous to Eldridge &Touil §004) as given in Table [3] The 
masses Mu(m,Z) and M core (m,Z) are dependent on metallicity. 
Therefore the lower and upper mass limits for these SN types are 
m SL,(i) = m SL,(i)(Z) and msu,(i) = ?«su,(i)(Z). The lower mass limit 
for Type II-P SNe is fixed at msL.iip(Z) = wsl.iip = 8 M . The 
absolute upper mass limit for the most massive stars is the cutoff 
mass defined by the IMF, thus msu,(i)(Z) = fnsu,® = m 2- 

The criterion for direct BH formation is based on the system 
used bv lHeger et all d2003l) and lEldridge & Toutl d2004l) . A direct 
BH forms when the final He core mass M COK (m,Z) exceeds the 
critical He core mass, M™^ = 15 M . Models using yields from 
IEIT08I will be referred to as 'EIT08M' . 

3.3.2. SN yields from rotating stellar models 

We take the yields for Mu(m,Z ), M fin (m,Z), F E (m,Z), and 
Yz(m, Z) from lGeorgy et al. ( 2009). Rotationally enhanced mass 
loss becomes very efficient with increasing m etallicity. The mass 
loss descri ption used in their models are from lMevnet & Maederl 
J20031 12001 with a rotational velocity of 300 km/s. Provided 
that a SN occurs even if a BH is formed, the resulting range for 
Type II SN (without further subdivision) is between 8-54 M at 
a metallicity Z = 0.004 Z . However, only stars between 8-25 
M Q form neutron stars. At Z the Type II mass range becomes 
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Table 2. List of all model parameters 



Parameters 


Value 


Unit 


Description 


Mm 


5 x 10 10 , 1 x 10", 5 x 10", 1.3 x 10 12 


M G 


Initial mass of the galaxy 


^ini 


1 x 10 3 


M Q yr" 1 


Star formation rate 


7 


lfr 6 


z Q 


Initial metallicity 


k 


1.5 




Power for the relation t//(f) oc Mism(0* 




800, 100, 




Swept-up ISM mass per SN 


M crit 

core 


15 


Ms 


Critical He core mass 




0.93 




SN dust destruction factor 


MsMBH 


3 x 10 9 , 5 x 10 9 


M G 


Mass of the SMBH 


's.MBH 


4 x 10 8 


yr 


Growth timscale for the SMBH 


'max 


10 9 


yr 


Maximum computed age of the galaxy 


Parameters 


Switch 




Description 


Y z , Y E , Y Q (for SN) 


EIT08, WW95. N06. Georgv et al. (2009) 




Possibilities for the SN yields 


Yz, Y E , Yq (for AGB) 


van den Hoek & Groenewegen (1997) 




Possibilities for the AGB yields 


<j>(m) 


Salpeter, mass-heavy, top-heavy, Larson 1, Larson 2 


Initial mass function 


SFR 


evolving / constant 




Additional switch for the SFR 


e AGB (m,Z) 


only one case considered 




Dust formation efficiency, AGB 




max / low 




SN dust formation efficiency 


fsN 


considered / not considered 




SN dust destruction 


BH/SN 


SN when BH / no SN when BH 




Possibility, if a SN occurs even a BH 








is formed or not 


SMBH 


considered / not considered 




Growth of SMBH 



narrower (8-25 M Q ) and stars above the upper limit will become 
WR stars and explode as Vole SNe. The range for neutron stars 
at this metallicity extends up to 35 M . 

To determine the lower msL,(i)(Z) and upper msu,(i)(Z) mass 
limits for the considered SNe types and remaining stars (i = IIP, 
II, Ib/c, R), we appled the same definitions as described for the 
EIT08M models. However, the direct BH cut is constrained by 
the mass of the remnant star M rem (m,Z) instead o f M c0le (m,Z). 
We fo llow the notation for the BH formation by Georg v et alj 
(120091) . where a BH forms when M lem (m,Z) > 2.7 M Q . Models 
where these yields and the above SN type division and BH for- 
mation are applied, are referred to as 'G09M'. 

3.3.3. Models with fixed SN mass range 

A treatment in which the boundaries of the mass range for 
all considered dust-forming supernovae and remaining stars are 
fixed throughout the evolutio n is quite common and has been 
used i n previous models (e.g.lDwek| [l 9 98; Morg an & Edmundsl 
l2003HDwek et al.ll2007r. IVahante et al.ll20(m 

We used the stellar yields for Y^(m,Z), Yz(m,Z) and certain 
elements Yg( m,Z) with Q = C,Q from nucleosyn thesis mod- 
els of either [Wooslev & Weaver! dl995L hereafter IWW95I) or 
iNomoto et alJ (120061) (herafter lNOfl) . All SNe are of Type II and 
considered to be within the mass interval otl,(H) = 8 M and 
m u,(ii) = 40. Yields for the final mass prior to explosion are not 
available, thus in all equations Fwind = 0. 

The remaining very massive stars between mL,(R) = 'Wu.(ii) = 
40 and the upper limit bi m] = m y, are assumed to turn into BHs. 
We apply the yields from EIT08 to these remaining stars in order 
to account for the gas return into the ISM from their mass loss 
phase according to Eqs.fTT1andfT2l Thus, the option of whether 
a SN explosion occurs or not is retained. 

The models with stellar yields from lWW95l wi ll be r eferred 
to as 'WW95M' and the models with yields from lN06l will be 
referred to as 'N06M'. 



3.3.4. AGB stars 

For high-mass AGB stars we use the stellar yields for the total 
ejected mass Y^(m,Z), the metal s Yjjm, Z) and certain elements 
Yo{m, Z) with Q = C,0 from Ivan den Hoek & Groenewegen! 
dl997l) . The mass interval is set to a fixed mass range of mL,(i) 
= 3 M and wu,© = 8 M Q , i = AGB. AGB stars with masses 
below 3 M R are not considered. 



3.4. AGB star and SN dust production efficiencies 

The values for the AGB star and SN dust production efficiencies 
are derived from observationally and theoretically determined 
dust yields. The dust production efficiency per stellar mass and 
metallicity, ej(m, Z), is calculated as 



£i(m,Z) — 



Y d (m,Z) 
Y z (m,Zy 



(21) 



where Y&(m, Z) is assumed to be the final amount of dust ejected 
into the ISM (i.e., produced and possibly processed through 
shock interactions). 

For AGB stars, we adopt ed for Yd(m,Z) the theo retically 
derived total dust yields from iFerrarotti & Gaill (2006) an d fo r 
Yz(m,Z) the yields from Ivan den H oek & Groenewegenl( ll997l) . 
We calculated the AGB star efficiency for different metallicities. 
For illustration, we show the metallicity-averaged efficiency, 
6agb('«X which is representative of the trend of the metallicity - 
dependent e\GB(m,Z) in Fig. Q] One notices that, on average, 
AGB stars between 3-4 M are the most efficient dust produc- 
ers. 

For the dust production efficiency for SNe we used two plau- 
sible limits: (i) an upper limit that was motivated by theoreti- 
cal SN dust formation models, and (ii) a lower limit based on 
dust masses inferred from observations of older SN remnants, 
such as Cas A, B0540-69.3 Crab nebula, and 1E 102.2-7219 
(e.g., iDouvion et al .] 1200 It iBouchet et all l2004t iGreen et al.1 
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Fig. 1. Dust-production efficiencies of massive stars. The dashed 
line represents the 'maximum' SN efficiency e mlLX (m) and the 
dashed dotted line the 'high' SN efficiency ehigh('«)- The solid 
line is the 'low' SN efficiency ei ow (m). The left dotted line rep- 
resents the metallicity-averaged AGB efficiency 6agb(w)- The 
vertical line marks the least massive (3 M ) AGB star consid- 
ered. 



2004 iKrause et al] 12004 iTemim et a J l2006t iRho et all [2008 



Sandstrom et al. 2008; Williams et al. 2008 UDunneet al.l l2009: 
Barlow et al. 2010). For the upper limit we applied two different 



comb inations of the du st and m etal yields, which ar e (i) Y^jm , Z) 
from iTodini & Ferraral d2001l) w ith YAm.Z) from IWW95I and 
(ii) YAm. Z) from lNozawa et all (l2003h with Y z (m, Z) from lN06l 
Averaging the obtained efficiencies over Z leads to the metal- 
licity independent 'maximum' SN dust production efficiency, 
e max (m). For the lower limit an average dust mass of 3 x 10~ 3 
M inferred from the SN remnants is applied to all SNe in the 
mass range 8^4-0 M Q . To calculate the SN dust p roduction effi- 
ciency, stellar yields from lWW95llN06l and lEIT08l are used. The 
efficiencies obtained with these yields are averaged, constituting 
the 'low' SN dust production efficiency, ei ow (m). 

Applying a dust destruction factor (see Sect. I2.21 i to e msLX (m) 
results in a reduced dust production efficiency. We refer to this 
efficiency as 'high' SN dust-production efficiency, ehigh('«) = £sn 
e max (m). The amount of dust (~ 2-6 x 10~ 2 M ) produced per 
SN with this efficiency is comparable to the higher dust masses 
derived in older SN remnants. 

We used a dust destruc tion coefficient, £sn = 0.93, following 
Bian chi & Schneider! (|2007). These SN efficiencies may possi- 
bly also arise from different scenarios of dust formation, such as 
grain growth and shock interactions in SN remnants. The aver- 
aged SN efficiencies e max (m) and ei ow (m), as well as the 'high' 
SN efficiency, 6high('«X are seen m Fig.Q] 

3.5. Dust destruction in the ISM 

Once a dust grain is injected into the interstellar medium, it is 
subject to either growth or to disruptive or destructive processes. 
We here focus on the destructive and disruptive ones due to su- 
pernova shocks. 

Disruptive processes are those that lead to fragmentation of 
large dust grains (radius > 1000 A) into smaller dust grains (ra- 
dius < 500 A). IJones et alj fl996) found that shattering due to 
grain-grain collisions dominates vaporization and therefore also 



determines the grain size redistribution, which is shifted towards 
smaller grains. 

Destructive processes return dust grains back into the gas 
phase. The main destruction of dust grains in the ISM comes 
from the sputtering caused by i nterstellar shock waves wit h 
shock velocities > 100 km s -1 dSeablll987t: IJones et ail [l 996). 
The destruction takes place because of high-velocity gas-grain 
impacts of smaller projectiles (radius < 100 A), such as ener- 
getic He + ions, onto dust grains. This results in the removal of 
dust species at or near the surface of the grains. 

A dust grain is exposed to thermal and nonthermal sputter- 
ing, as well as to v aporization during the passage of the shock. 
IJones et all l[l996) also found that graphite grains are mainly 
destroyed by thermal sputtering. Silicates are equally affected 
by thermal and non-thermal sputtering in high-velocity shocks 
with velocities v, >150 km s , while vaporization is negligi- 
ble. These processes determine the lifetime of dust grains in the 
ISM. The timescale Tdi(f) of the dust grains against destruction, 
or simply the lifetime of the dust gra ins, and the dust destruction 
rate in the ISM £ D ,i SM (0 are given dMcKedl 19891 iDwekll 19981: 
iDwek et alJl200 7) through 



Tdl(0 

£d,ism(0 



M lSM (t) 
M cX R SN (t) ' 
MS) 

Tdl(0 ' 



(22) 



where Mism(0 is the mass of the ISM and /?sn(0 the supernova 
rate of all supernovae causing the destruction. The mass M c i is 
the mass of the ISM, which is completely cleared of dust through 
one single supernova remnant. These two equations combined 
lead to an expression for £ , d,ism(0 in the form 



£d,ism(0 = m(t)M cl R SN (t), 



(23) 



which shows that besides /?sn(0 and M c \, the dust destruction 
rate in the ISM also depends on the dust-to-gas ratio rjd{t). 

An expression for M c \ whi c h is d ependent on the shock ve- 
locity is given bv IDwek et ail (120071) . For a homogenous ISM 
and under the assumption t hat si licon and carbon grains are 
equally mixed, iDwek etalJ d2007l) obtains M cl = 1100-1300 
M . However, the ISM is inhomogeneous, charact erized by cold, 
warm, and hot phases with different densities (e.g.. lMcKedl 989, 
and references therein). The density contrast between the cold 
and warm phases and the hot phase can be relatively large. 
Shocks traveling through these phases are found to be very inef- 
ficient in destroying dust dJone s 2004). The destruction process 
is solely effective in the warm (T > 100K) phase of the interstel- 
lar medium, while SN shocks propagating through either a hot 
ISM with low density or cold clo uds (atomic a nd molecular) do 
not destroy dust effectively (e.g., McKea ll989l) . 

Another important parameter is the injecti on timescale o f 
stellar yields and dust into the ISM. Following [McK ee (1989), 
we estimate the injection timescale of the dust from stellar 
sources as 



Tin(0 



M d (0 



£d,SN(0 + £d,AGB(0 



(24) 



In GAH1 1 we made a rough estimate of the minimum averaged 
dust injection rate from SNe and AGBs of 0.5 M yr~' based on 
data for high-z QSOs. Using this dust injection rate we obtained 
an average injection time T m (t) = 4 x 10 8 yr. 

For comparison we estimate the lifetime Tdi(f) of the dust 
grains. The mass of the ISM is assumed to be Mism = 2 x 10 10 
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M Q . The SN rates are calculated for a constant SFR of 500 M 
yr -1 and for a Larson 2 IMF. This results in a timescale of Tdi(f) 
~ 1 .3 x 10 7 yr for M d = 100 M . A shorter timescale Tdi(0 ~ 1 .6 
X 10 6 yr is obtained for M d = 800 M . When using a Salpeter 
IMF for the SN rates, the timescales are usually longer (Tdi(f) ~ 
6.6 x 10 7 yr, for M d = 100 M ). 

We note that this is a rather rough estimate. Typically Ti n (t) 
and Tdi(?) are strongly dependent on the IMF, Mism(0 and the 
SNe dust production efficiency e;(m, Z). Hence, the values will 
deviate from the above estimated average during evolution. 

However, this example demonstrates that the dust injection 
timescale can be longer than the lifetimes of the dust grains. The 
difference between these timescales is influenced by the value of 
M c b Pertaining to the formation of high dust masses in galaxies, 
this has significant consequences. An injection timescale Ti n (f) 
longer than the destruction timescale Tdi(f) does not allow the 
build up of high dust masses in the galaxy. This implies that 
the dust injection rate must be higher than the dust destruction 
rate. A lowering of the dust destruction rate could be achieved if 
dust grains are shielded from destruction. Alternatively, a rapid 
dust grain growth in the ISM might be an option, if SNe and 
AGBs cannot generate the necessary high dust injection rate. 
Grain growth, however, is not incorporated into our model and 
remains to be investigated. 

Furthermore, we assume that the starburst occurs in an ini- 
tially dust free galaxy. Consequently dust produced by the first 
generations of SNe might be ejected into the ISM unhindered. 
The epoch at which the first SN shocks are able to sweep up 
ISM gas mixed with dust remains elusive. In view of these con- 
siderations, we stress that dust destruction in the ISM is uncer- 
tain, particularly when considering galaxies with conditions as 
described in Sect. 13. Il and l3.2l 

Despite our simple assumption of dust being immediately 
homogeneously mixed with the gas in the ISM, we account for 
the uncertainty of the lifetime of dust grains against destruc- 
tion by using M c \ as a parameter. This has also been done by 
iDwek et aH (120071) . The considered cases are for M c i = 800 M 
as the highest destruction, M c \ =100 M for modest destruction, 
and M c i = for no dust destruction. 



4. Results 

In this section we present the results of several models calcu- 
lated within the first Gyr after starburst. We focus on the total 
dust mass in a galaxy. At redshift z > 6, the time to build up high 
dust masses of > 10 8 M is limited to 400-500 Myr. Models 
able to exceed an amount of 10 8 M of dust within this time are 
therefore of particular interest. General evolutionary tendencies 
of certain quantities, such as dust injection rates, SFR, metallic- 
ity, or the amount of gas, are also discussed. In this paper we 
present results that assume the sa me value for the i nitial SFR in 
all models as given in Table |2 In lGall et alJ d201 ll) . we investi- 
gated models for different initial SFRs and discuss their applica- 
bility to particular quasars at z > 6. 

4. 1. Evolution of dust, dust-to-gas, and dust-to-metal 
ratios 

We study models with all considered initial gas masses M m \. The 
models discussed are calculated for the case that no SN occurs 
when a BH is formed (see Sect. l2.2.lT ). Models including SMBH 
formation are deferred to Sect. I4.5I 



4.1 .1 . Models with SN type differentiation 

EIT08M with Mi ni = 5x10" M : In Fig.PJthe results are pre- 
sented for the mass of dust Md(f) in a galaxy with M- m \ = 5 x 10 11 
M . Calculations were performed for a 'maximum' (top row), a 
'high' (middle row), and a 'low' (bottom row) SN dust produc- 
tion efficiency e;(m). Dust destruction in the ISM is included for 
three values of M c i (see Sect. 13.51 ). This model is considered as a 
reference model. We therefore describe it in detail. 

Dust masses Md(f) obtained with the 'maximum' SN dust 
production efficiency e max (m) exceed 10 8 M of dust in all cases 
of M c i and for all IMFs. In fact, 10 8 M of M d (f) is reached 
within the first few Myr. This is insensitive to both the IMF and 
the amount of M c i. The latter however determines the suppres- 
sion of M d (f) as the system evolves. For M d = 800 M (top left), 
a maximum amount of ~ 4 x 10 8 M is reached for a Larson 
2 IMF and sustained until the end of the computation. Only the 
IMFs favouring lower masses exhibit a shallow decline in Md(?)- 
Dust destruction with M c i =100 M (top middle) leads to nearly 
constant dust masses retainable over 1 Gyr of evolution. The 
amount of dust for most IMFs is a few times 10 9 M . When 
assuming no dust destruction (top right), a broader spread of the 
dust masses for the different IMFs develops. The amount of dust 
increases with time regardless of the IMF and yields up to 10 9 ~ 10 
M are reached. 

In the case of 'high' SN dust production efficiency, e n igh('«X 
similar trends are featured. However the amount of dust is lower. 
Without dust destruction (M c i = 0), dust masses up to a few 
times 10 8 M are attained for all IMFs (Fig. [2] middle right). 
The timescale to exceed 10 8 M of dust ranges from 80 Myr 
(Larson 2 IMF) up to 400 Myr (Salpeter IMF). Taking dust de- 
struction with only a modest amount of M c i = 100 M into ac- 
count decreases Mi(f) substantially (Fig. [2] middle middle). The 
IMFs favouring high masses are most affected, while the reduc- 
tion in Md(f) for a Salpeter IMF is small. Except for the Salpeter 
IMF, all remaining IMFs lead to more than 10 8 M of dust. The 
highest amount of dust is reached with either a Larson 2 or a 
mass-heavy IMF, and is ~ 2 x 10 8 M at an epoch of ~ 400 Myr. 
Considering a very high destruction with M c \ = 800 M results 
in a strong reduction of Md(f), so only ~ 2-3 x 10 7 M of dust 
remains throughout the evolution (Fig. [2] middle left). 

The results for Md(f) with a 'low' SN dust production ef- 
ficiency ei ow (m) are dominated by AGB dust production (see 
also GAH11). Without dust destruction, dust masses up to 1— 
2 x 10 8 M is reached after ~ 300-500 Myr for most of the 
IMFs. Applying a modest dust destruction of M c i = 100 M re- 
duces Mi(f) analogously to the higher SN efficiencies e max (m) 
or ehigh('w)- Dust masses for either the top-heavy or the Larson 
2 IMF are efficiently decreased and comparable to Md(?) for a 
Salpeter IMF. Increasing M c i to 800 M leads to a stronger re- 
duction in Md(f) for these IMFs, resulting in lower dust masses 
than for a Salpeter IMF. 

A considerable difference in the progression of Md(f) for the 
'low' SN efficiency ei ow (m) compared to the higher SN efficien- 
cies is encountered during the first ~ 200 Myr. While a fast rise 
of Md(f) is identifiable for either e max (m) or ehigh^X f° r ei ovl (m) 
the dust mass remains between 10 6 ~ 7 M during this epoch. This 
is caused primarily by high-mass AGB stars with short lifetimes 
(> 4-5 M ). A further increase in the dust mass due to the de- 
layed AGB dust injection from the less massive but more effi- 
cient AGB stars results in a second dust bump at an epoch of ~ 
200-300 Myr. High dust masses at early epochs (100-200 Myr) 
are not possible with ei ow (m). 
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Fig. 2. Evolution of the total dust mass and dust destruction rates for EIT08M. The initial gas mass of the galaxy M m \ = 5 x 10 11 
M . Calculations are performed for a 'maximum' (top row), a 'high' (middle row), and a 'low' (bottom row) SN dust production 
efficiency £i(m). Dust destruction is taken into account for M c \ = 800 M (left column), 100 M (middle column), and M (right 
column). The bottom group of curves in the light blue area represents the total dust destruction rate E£,(t), which is based on Eq.[T6l 
The upper group of curves in the green or yellow zones displays the evolution of the total amount of dust Md(f) in the galaxy. The 
yellow area marks dust masses exceeding 10 8 M of dust. The grey horizontal dashed line marks the limit of 10 8 M of dust. The 
grey vertical dashed lines indicate epochs at 100, 170, 400, and 500 Myr after the onset of starburst. The black solid, green dotted, 
cyan dashed, magenta dashed-dotted, and blue dashed-dot-dotted lines represent the Salpeter, mass-heavy, top-heavy, Larson 1 and 
Larson 2 IMF, respectively. 



The curves at the bottom of the light-blue area in Fig. [2] 
represent the dust destruction rate Eo(f), F° r M c \ = and in- 
dependent of e\(m) the dust destruction rate reflects the amount 
of dust incorporated into stars per unit time and is calculated as 
E D (t) = j] d (t) iff(t). For M c i > E D (t) is given by Eq. [16] and 
additionally consists of the destroyed dust in the ISM, #d,ism(*). 

We find that Eo(t) increases faster and earlier with increas- 
ing M c \. Consequently, dust destruction rates comparable to SN 
injection rates are reached at earlier epochs. Given Eq.|23]for the 
amount of dust destroyed through SN shocks, the high SN rates 
RsN(t) at the beginning of evolution in conjunction with M c \ lead 



to a high base value of £ , d,ism(0- This suppresses the rise in the 
amount of dust at early epochs. Additionally, SN rates are higher 
for the IMFs biased towards more massive stars affecting these 
IMFs most. Later in the evolution, the dust destruction rates re- 
main almost constant in most cases, with a marginal increase 
for M c i = or decline for M c i > 0. Dust destruction and injection 
appear to be balanced, resulting in the flat development of M^(t). 

In Fig. [3] we present the results for the evolution of the dust- 
to-gas mass ratio r]d(i) and the dust-to-metal mass ratio rjzd(t)- 
Commonly the curves for T)&(t) for all IMFs are slowly increas- 
ing with time except for M c i = 800 M for which they remain ap- 
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Fig. 3. Evolution of the dust-to-gas and dust-to-metal mass ratios for E1T08M. The initial gas mass of the galaxy M; n i = 5 x 10 11 
M . Calculations are performed for a 'maximum' (top row), a 'high' (middle row), and a 'low' (bottom row) SN dust production 
efficiency £i(m). Dust destruction is taken into account for M c \ = 800 M (left column), 100 M (middle column), and M (right 
column). The upper group of curves signifies the dust-to-metal mass ratio //zd(0- The lower group of curves represents the dust- 
to-gas mass ratio rj&it). The gas-to-dust ratio displayed in the figures is calculated for a Larson 2 IMF at an epoch of 400 Myr. 
The grey vertical dashed lines indicate epochs at 100, 400, and 500 Myr after the onset of starburst. The black solid, green dotted, 
cyan dashed, magenta dashed-dotted, and blue dashed-dot-dotted lines represent the Salpeter, mass-heavy, top-heavy, Larson 1 and 
Larson 2 IMF, respectively. 



proximately constant. Interestingly, in this case and for e max (m), 
the dust-to-gas ratio for all IMFs sustains a constant value of 
~ 10~ 3 over the whole evolution. Without dust destruction (top 
right), 77d(0 increases up to values above ~ 10~ 2 for most of 
the IMFs. For either a 'high' or 'low' SN efficiency, the values 
are below ~ 10~ 3 . Only for ehi g h(w) and M c \ = (middle right) 
is ?7d(f) further increased. These overall trends show that rj&(f) 
can be very low for massive galaxies even if the galaxy appears 
dusty. 

The dust-to-metal mass ratio rjzdit) is significantly lower for 
the Larson 2 and top-heavy IMF than for the other IMFs. This 
feature is exhibited in all calculated models. The degree of this 



separation between different IMFs depends primarily on the SN 
dust production efficiency and secondly on the destruction rate 
in the ISM, and is larger for ei ow (m). This reflects the decrease in 
the SN dust production efficiency towards the higher mass end 
of the SN mass interval. Additionally, Type Vole SNe produce 
no dust, but they do inject metals into the ISM at higher rates 
for the top-heavy IMFs. For M c \ = and ehigh('«X the amount 
of metals bound in dust grains is between ~ 1.5-3 %. For M c \ 
= 100 M the fraction of metals bound in dust grains for the 
two IMFs favouring high masses is below 1 % and decreases 
further with increasing M c \. Generally, for a 'low' SN efficiency, 
rjzd(t) remains below 10~ 2 for all IMFs. For the 'maximum' SN 
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Fig. 4. Evolution of the total dust mass for EIT08M. Calculations for galaxies with initial gas masses of M m \ =5 x 10 M G are 
presented in the first and second columns and for M m \ = 1 x 10 11 M in the third and fourth columns. Total dust masses M&{t) 
are shown for a 'maximum' (top row), a 'high' (middle row), and a 'low' (bottom row) SN dust production efficiency €\(m). Dust 
destruction is taken into account for M c \ = 100 M (first and third columns) and M (second and fourth columns). Maximum dust 
masses exceeding 10 8 M are indicated as yellow shaded zones. The grey horizontal dashed line marks the limit of 10 8 M of dust. 
The grey vertical dashed lines indicate epochs at 400, and 500 Myr after the onset of starburst. The black solid, green dotted, and 
blue dashed-dot-dotted lines represent the Salpeter, mass-heavy and Larson 2 IMF, respectively. 



efficiency e mdx (m) the dust-to-metal ratio, r]zd(t), is about a factor 
of 10 higher than for ehigh('w). 

EIT08M with M ini = 1 x 10 11 M : In Fig.g](two right columns) 
the results are presented for Md(f) in cases of either no or a 
modest destruction and all three SN dust formation efficiencies 
emax(ra), en S h(m), e\ ow (m). 

The top row of Fig. [4] depicts the evolution of the amount of 
dust Md(f) for € mm (m). It is evident that high dust masses beyond 
10 8 M are obtained for both cases of M c i. The maximum value 
for Md(f) is already reached shortly after the onset of the star- 
burst. Thereafter Md(f) follows a negative slope, which is steep- 
est for the IMFs biased towards low-mass stars. This decline is 
generally observable in all cases of M c \ and e;(m). The cause of 
this is that dust is incorporated into stars on higher rates than 
is replenished by stellar sources. This dust decrease is amplified 
by dust destruction in the ISM as is demonstrated for M c \ = 100 
M . For e n igh(?w) and M d = 0, a high total dust mass > 10 8 M 



can be achieved with a Larson 2 IMF after ~ 200 Myrs (Fig. [4] 
middle row, right column). For all remaining cases of either M c \ 
or for ei ow (m) and regardless of the IMF, M&(f) stays below 10 8 
M . The tendencies for the various IMFs resemble those identi- 
fied for M ini = 5 x 10 11 Mo. 

The curves of the dust-to-gas ratio rji(t) and the dust-to-metal 
ratio ?7zd(0 are also similar to the system with Mj n ; = 5 x 10 
M , although the amount is higher (see Fig.|5]two left columns). 
For M c i = and 'maximum' SN efficiency, rja{t) is roughly inde- 
pendent of the IMF, while exhibiting a wider range for ei ow (m). 
The fraction of metals bound in dust grains is < 4 % (M c i = 0) 
for IMFs favouring low-mass stars and ehigh(^); for e m:ix (m) it is 
roughly 20-30%. 

EIT08M with Mini = 5 x 10 10 M : The evolution of M d (f) 
within the first Gyr exhibits similar, but more strongly pro- 
nounced, tendencies for the various IMFs as the galaxy system 
with Mi n i = 1 x 10 11 M . Results are presented in Fig.|4](two left 
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Fig. 5. Evolution of the dust-to-gas and dust-to-metal mass ratios for EIT08M. Calculations for galaxies with initial gas masses of 
Mi n i = 5 x 10 10 M are presented in the first and second columns and for M m \ = 1 x 10 1 1 M in the third and fourth columns. Results 
are presented for a 'maximum' (top row), a 'high' (middle row), and a 'low' (bottom row) SN dust production efficiency e;(m). Dust 
destruction is taken into account for M c \ =100 M (first and third columns) and M (second and fourth columns). The upper group 
of thick curves signifies the dust-to-metal mass ratio rjzd(t). The lower group of thin curves represents the dust-to-gas mass ratio 
T]i{t). The gas-to-dust ratio displayed in the figures is calculated for a Larson 2 IMF at an epoch of 400 Myr. The grey vertical dashed 
lines indicate epochs at 400, and 500 Myr after the onset of starburst. The black solid, green dotted, and blue dashed-dot-dotted 
lines represent the Salpeter, mass-heavy and Larson 2 IMF, respectively. 



columns). Even with a 'maximum' SN efficiency e max (m), large 
dust masses of a few times 10 8 M cannot be reached or sus- 
tained when dust destruction in the ISM is considered. In Fig. [5] 
(two left columns) it is evident that the dust-to-gas ratios 7/d(f) 
and the dust-to-metal ratios 77zd(0 result in higher values than 
for the higher mass galaxy systems (see Figs. [3] [5] two right 
columns). The tendencies are the same in general. 



G09M: In Fig. |6]the results are shown for models with a rota- 
tionally enhanced mass-loss prescription. The galaxy under con- 
sideration has Mi n i = 5 x 10 11 M . The outcome is similar to the 
reference EIT08M (Fig.f2]i, although the obtained dust yields are 
increased for the IMFs favouring higher masses. For a Larson 2 
IMF, M d (f) exceeds 10 9 M after ~ 650 Myr. Interestingly, even 
though the dust yields are high for no dust destruction in the 
ISM, for M c i = 800 M the dust mass also stays below 10 s M . 
In case of a 'low' SN dust production efficiency, 10 8 M cannot 



be reached within the first 400 Myr. However, for the IMFs bi- 
ased towards high masses an amount of dust > 10 8 M is attained 
later in the evolution. 



4.1.2. Models with fixed SNe mass range 

In Fig.|7]the results of M d (f) are shown for WW95M (top row) 
and N06M (bottom row) for a galaxy with M in i = 5 x 10" 
M and 'high' SN efficiency e n igh('«)- The amount of dust in 
N06M early in the evolution is greater than for WW95M for 
the IMFs favouring massive stars, but flattens later, leading to 
slightly lower dust masses than achieved with WW95M. The 
dust masses achieved with either the 'maximum' or the 'low' 
SN dust production efficiencies are nearly identical to EIT08M 
and are therefore not shown. Generally, the evolution of these 
models is similar to EIT08M. 
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Fig. 6. Evolution of the total dust mass for G09M, presented for 
a galaxy with M iai = 5 x 10 11 M . Total dust masses M d (f) are 
shown for a 'high' (top row) and a 'low' (bottom row) SN dust 
production efficiency e;(m). Dust destruction is taken into ac- 
count for M c i = 800 M (left column) and M G (right column). 
Maximum dust masses exceeding 10 8 M are indicated as yel- 
low shaded zones. The grey horizontal dashed line marks the 
limit of 10 8 M of dust. The grey vertical dashed lines indicate 
epochs at 400 and 500 Myr after the onset of starburst. The black 
solid, green dotted, cyan dashed, magenta dashed-dotted, and 
blue dashed-dot-dotted lines represent the Salpeter, mass-heavy, 
top-heavy, Larson 1, and Larson 2 IMF, respectively. 
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Fig. 7. Evolution of the total dust mass for WW95M (top row) 
and N06M (bottom row). Both models are presented for a galaxy 
with M mi = 5 x 10 11 M . Total dust masses M a (t) are shown for 
'high' SN dust production efficiency ehigh('w)- Dust destruction 
is taken into account for M c i = 100 M (left column) and M 
(right column). Maximum dust masses exceeding 10 8 M are 
indicated as yellow shaded zones. The grey horizontal dashed 
line marks the limit of 10 8 M of dust. The grey vertical dashed 
lines indicate epochs at 400 and 500 Myr after the onset of 
starburst. The black solid, green dotted, cyan dashed, magenta 
dashed-dotted, and blue dashed-dot-dotted lines represent the 
Salpeter, mass-heavy, top-heavy, Larson 1, and Larson 2 IMF, 
respectively. 




4.1 .3. Very high M ini and the case of constant SFR 

In Fig. [8] the results are displayed for a case with constant SFR 
if/(t) = 1 x 10 3 M yr and a galaxy with a high initial gas mass 
of Mi n i = 1.3 x 10 12 M . Such a high M; n j is necessary when 
assuming a constant SFR, since in the lower mass galaxies the 
available mass for star formation gets exhausted before an age 
of 1 Gyr is reached. 

We find that a build up of high dust masses also for such 
massive galaxies is greatly suppressed with a high amount of 
destruction. With none of the IMFs, 10 8 M of dust is attained, 
apart from the models with e max . High dust masses with either 
thigh or eiow are only possible in case of modest or no dust de- 
struction. For the latter the results are shown in the right column 
of Fig. [8] Interestingly, also for 'low' SN dust production effi- 
ciency, dust masses > 10 8 M are possible after ~ 300^1-00 Myr. 

The results are very similar to an EIT08M with such high 
Mi n i and an evolving SFR. A slight difference appears for M^it) 
and IMFs that favour low-mass stars. These models remain close 
to constant after 500-600 Myr, while the models for constant 
SFR slowly decline. The similarity between these models can be 
explained by the very high mass of the galaxy. In neither model 
does the ISM mass get significantly reduced. 



4.2. Evolution of dust production rates, SFR and 
metallicity 

In Fig. [9] we present the evolution of quantities such as the total 
dust injection rates E a ^(t), £d,AGB(0i the AGB and SNe rates 
^agb(?X ^sn(?X the SFR and the metallicity Z(f) for a range of 
initial galaxy gas masses. The rates and the metallicities are in- 
dependent of the assumed dust formation efficiency and can be 
discussed for each initial mass Mj n ;. The results for these quan- 
tities, which come from EIT08M, WW95M, and N06M with 
identical parameter setting, are very similar. Consequently these 
quantities are discussed based on EIT08M. The upper row of 
Fig.|9]shows the dust injection rates of AGB stars £d,AGB(0 and 
SNe £ , < i,sn(0- For supernovae £ , ( i,sn(0 is highest for a Larson 2 
IMF and lowest for a Salpeter IMF. 

The SN dust production rates are approximately one order 
of magnitude lower for ei ow (m) than for ehigh('«)- Similar differ- 
ence is the case between £hi g h(m) and e mdx {m). In all cases a clear 
separation of the values of £ , ( i,sn(0 amongst the various IMFs is 
pronounced throughout the evolution. 

The dust injection rates for AGB stars, £d,AGB(fX are con- 
siderably influenced by the long lifetimes. After ~ 200 Myr the 
AGB dust production rates M^agbW are comparable to the dust 
injection rates for SNe with ehigh('«)- This is caused by the higher 
dust production efficiency for AGB stars between ~ 3-4 M , 
leading to larger amounts of dust produced than for AGB stars 
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Fig. 8. Evolution of the total dust mass for WW95M with con- 
stant SFR. WW95M is presented for a galaxy with M ini = 1.3 x 
10 12 M . Total dust masses Md(f) are shown for 'high' (top row) 
and a 'low' (bottom row) SN dust production efficiency e;(m). 
Dust destruction is taken into account for M c i = 800 M Q (left col- 
umn) and M (right column). Maximum dust masses exceed- 
ing 10 8 M are indicated as yellow shaded zones. The grey hor- 
izontal dashed line marks the limit of 10 8 M of dust. The grey 
vertical dashed lines indicate epochs at 400 and 500 Myr after 
the onset of starburst. The black solid, green dotted, cyan dashed, 
magenta dashed-dotted, and blue dashed-dot-dotted lines repre- 
sent the Saipeter, mass-heavy, top-heavy, Larson 1, and Larson 
2 IMF, respectively. 



in the mass range of 4-8 M (see Sect. 13.41 ). The small variations 
in £d,AGB(f) are mainly due to alteration of the stellar yields and 
dust formation efficiencies at different metallicities. 

The SNe and AGB dust injection rates decline as a conse- 
quence of decreasing SNe and AGB rates, as shown in Fig. [9] 
(second row). Both rates exhibit a faster decline for the lower 
massive galaxies and for IMFs biased towards the intermediate 
and low mass stars. 

This behaviour is determined by the slope of the SFR, which 
depends on the initial mass of the galaxy (see Fig. |9]third row). 
The decline in the SFR for systems with M- mi > 3 x 10 11 M 
is shallower. For these galaxies, a high SFR of a few hundred 
M yr can be sustained over at least lGyr of evolution. For 
the lower mass galaxies, the SFR declines faster and exhibits a 
strong dependence on the IMFs. At a time of 400 Myr the SFR 
for a galaxy with M ini = 1 x 10 11 M is between 60 (Saipeter) 
to 500 (Larson 2) M yr" 1 . When M mi = 5 x 10 10 M , if/(f) = 
20^1-00 M yr~' . After ~ 1 Gyr the difference between the SFRs 
obtained with a Larson 2 IMF and a Saipeter IMF is more than 
an order of magnitude. 

In Fig. [9] (bottom row) it is seen that the metallicity Z(f) in 
the two lower mass galaxies rises quickly within the first 100- 
200 Myr up to values of more than 5 Z for either a top-heavy 
or a Larson 2 IMF. Thereafter it remains rather constant. For 
the system with M; n ; = 5 x 10 11 M , the metallicity increases 



more slowly than in the lower mass galaxies. At an epoch of 400 
Myr, a metallicity of 2-3 Z can be observed for the top-heavy 
and Larson 2 IMFs, while the Saipeter IMF achieves only a bit 
less than half solar. In the most massive galaxies, the metallicity 
exceeds Z only in case of the top-heavy IMFs before an age 
of 400-500 Myr, but remains below Z for the IMFs favouring 
lower mass stars. 

4.3. Evolution of gas, metals and stellar masses 

The evolution of quantities, such as the gas mass, mass of metals, 
and stellar masses is discussed based on EIT08M. The Fig. [10] 
illustrates the evolution of the H + He gas mass Mc(f). As a 
consequence of the scaling of the SFR with the mass of the ISM, 
M IS m(0 = M G {t) + M z (t), the progression of M G {t) is identical 
to the SFR described in Sect. 072] 

The most massive galaxy with M m i = 1.3 x 10 12 M has a 
residual gas mass Mc(t) of ~ 10 12 M after 1 Gyr. This explains 
the very low dust-to-gas ratio for such a system even though 
very high dust masses can be reached. The galaxy with Mi„i = 
5 x 10 11 M gets more exhausted, but also retains a gas mass 
M G (?) of 3^ x 10 1 1 M after ~ 400-500 Myr. 

In analogy with the curves of the SFR, the curves of the gas 
mass for the system with M m \ = 1 x 10 11 M feature a pro- 
nounced separation between the IMFs. A Larson 2 IMF shows 
a flat evolution leading to a slower depletion of the gas than for 
the Saipeter IMF and after 400-500 Myr ~ 5 x 10 10 M of gas 
is still present. The lowest mass system (M m \ = 5 x 10 10 M ) 
exhibits the strongest dependence on the IMF. After 400 Myr, 
M c (t) is ~ 2 x 10 10 M for a Larson 2 IMF. The difference be- 
tween a Saipeter IMF and a Larson 2 IMF after one Gyr is an 
order of magnitude. The residual gas mass with a Saipeter IMF 
is a few times 10 8 M at this time. 

The steepness of the decline in the progression of Mc(t) is 
influenced by the stellar feedback. For either a Larson 2 or a top- 
heavy IMF, more massive stars are formed. Such stars live short 
lives and release most of their mass back into the ISM, through 
either stellar winds or in explosive events. The IMFs favouring 
lower masses, however, lock most of the gas used for star forma- 
tion into intermediate-to-low-mass stars. These stars are formed 
at higher rates than massive stars. The low-mass stars live a long 
time and also do not inject much material back into the ISM dur- 
ing our considered time span of 1 Gyr. As a consequence, the 
material available for star formation gets more rapidly depleted 
for the IMFs favouring low-mass stars than for the top-heavy 
IMFs. Hence, the evolution of the gas mass and the SFR results 
in a steeper decline. 

This contrasts with the stellar masses M„(f) (Fig.[TT]top row) 
obtained from the relation M»(f) = M ini - M IS m(0- For the IMFs 
favouring lower mass stars, M»(f) rises steeply and approaches 
the initial mass of this galaxy in the case of the lower massive 
galaxy systems. The slope of M»(f) for the top-heavy IMFs is 
shallower, and lower stellar masses M,(/) are achieved. 

The Fig. QT| ( bottom row) depicts the mass of the ejected 
heavy elements M z (f). For systems with M ini > 1 x 10 11 M , the 
metal enrichment in the ISM is considerable. The mass of the 
heavy elements increases up to a few times 10 10 M . 

In the lower mass galaxies M m \ < 1 x 10 11 M , the amount 
of metals Mz(t) attained reaches a maximum within the first 200 
Myr whereafter Mz(f) declines. This is caused by astration. The 
maximum mass of metals obtained with a Saipeter IMF is only 
a few times 10 8 M . For the lowest mass galaxy the amount 
of metals after ~ 600 Myr for a Saipeter IMF is less than 10 8 
M . This implies that dust masses Md(f) in excess of 10 8 M for 
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Fig. 9. Evolution of total dust injection rates for AGB stars and SNe, the AGB and SNe rates, the SFR, and the metallicity. Results 
are shown for M ini = 1.3 x 10 12 M Q (first column), 5 x 10 11 M (second column), 1 x 10 11 M (third column), and 5 x 10 10 
M (fourth column). The solid, dotted, dashed, dashed-dotted, and dashed-dot-dotted lines represent the Salpeter, mass-heavy, top- 
heavy, Larson 1, and Larson 2 IMF, respectively. First row: Total SN dust injection rates £<i,sn(0 f° r a 'low' (black lines), 'high' 
(blue lines), and 'maximum' (green lines) SN efficiency e;(m), and AGB dust injection rates £d,AGB(0 (red lines). SN dust injection 
rates are only shown for a Salpeter (solid) and a Larson 2 (dashed-dot-dotted) IMF. Second row: SNe rates /?sn(0 (black lines) and 
AGB rates /?agb(0 (red lines). SN rates are only shown for a Salpeter (solid), mass-heavy (dotted), and a Larson 2 (dashed-dot- 
dotted) IMF. Third row: Evolution of the SFR. The blue area marks the region between a SFR of 100-1000 M Q yr~'. Fourth row: 
Evolution of the metallicity Z(t). The green regions mark a metallicity between 1-2 Z (light green), 2-3 Z (dark green), 3>-A Z 
(grass green). The grey vertical dashed lines indicate epochs at 400 and 500 Myr after the onset of starburst. 



this IMF are unfeasible, even if dust-grain growth in the ISM is 
invoked. 



4.4. Comparison to other dust evolution models for high-z 
galaxies 

The models with a 'maximum' SN dust production efficiency for 
a system with M ini of 5 x 10 10 M (Fig.|4]two left columns, first 



row) can be compared to models of iDwek et al. (2007). For this 
efficiency the AGB dust production is neglig ible and the SN dus t 
yields are similar to the dust yields used bv lDwek et al.l J2007). 
Our results for a top-heavy IMF and M c \ = at an epoch of 400 
Myr agrees with their results, while they disagree for a Salpeter 
IMF and for cases with M d = 100 M Q . 

The origin of this disagreement can be traced to the neglect 
of the lifetime-dependent mass injection of stars and a differ- 
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Fig. 10. Evolution of the gas mass M G {t) based on EIT08M. 
Results are shown for M ini = 1.3 x 10 12 M (black lines), 5 x 
10 11 M (red lines), 1 x 10 11 M Q ( green lines), and 5 x 10 10 M G 
(magenta lines). The grey vertical dashed lines indicate epochs 
at 400 and 500 Myr after the onset of starburst. The solid, dotted, 
dashed, dashed-dotted, and dashed-dot-dotted lines represent the 
Salpeter, mass-heavy, top-heavy, Larson 1, and Larson 2 IMF, 
respectively. 



4.5. Results for models including a SMBH 

We have investigated whether including the formation of the 
SMBH leads to differences in the evolution of either the total 
amount of dust or the physical properties of a galaxy. Based on 
the observed SMBH masses > 10 9 M r for high-redshift QSOs 
(e.g.. IWillott et all 12003b IVestergaardlT2004fc I Jiang et alJ 120061: 
IWang et al.1 1201 Oh . we studied two cases for the SMBH mass, 
M SM bh = 3 x 10 9 M an d M S mbh = 5 x 10 9 M Q . According to 
Kawakatu &Wadal (l2009h . the final SMBH mass takes up ap- 
proximately 1-10 % of the supply mass Msmbhsup- This implies 
that to grow a SMBH of 3(5) x 10 9 M the minimum required 
supply mass ranges between 3(5) x 10 10 to 3(5) x 10 11 M G . In 
view of our assumption that Msmbhsup = Mjm, we included the 
SMBH formation in galaxies with M ini = 5 x 10 10 M and M ini 
= 1 x 10 11 M . 

In Fig. [12] we show the results of Md(f) for a Salpeter and 
Larson 1 IMF for an EIT08M and M ini = 5 x 10 1() M where the 
SMBH formation has been included. These results are compared 
to those of a model with the same M m i, but without a SMBH. 

Taking the SMBH into account leads to a steeper decline in 
Md(f) and therefore to a lower amount of dust than for the model 
without a SMBH after ~ 100 Myr. All models are similar within 
the first ~ 100 Myr. For M SM bh = 5 x 10 9 M q the difference is at 
most ~ 50 % for a Salpeter IMF and ~ 30 % for a Larson 1 IMF 
at an epoch of 400 Myr. At the same epoch but for Msmbh = 3 
x 10 9 Mo, the dust mass for a Salpeter IMF is reduced by only 
about 30 %, and 20 % for the Larson 1 IMF. In cases of IMFs 
biased towards higher masses and in the more massive galaxy 
systems, the formation of the SMBH, as introduced here, does 
not noticeably effect the progression of Md(f) and the physical 
properties of a galaxy. 



en t treatment of the mass recycled into the ISM in the models 
of lDweketai](l2007h . The latter is approximated as 0.5if/(t) in 
their models, which leads to an IMF-independent evolution of 
the ISM mass. In the case of a Salpeter IMF, this simple ap- 
proximation has consequences for the evolution of all physical 
properties. It implicitly presupposes that even stars, which are 
formed between the lower mass limit of the IMF and the lower 
SN mass limit, immediately return half of their stellar mass back 
into the ISM. However, stars between 3-8 M eject their ele- 
ments up to a few 100 Myr delayed and stars < 3 M do not 
eject a significant amount of elements within the first Gyr. For a 
Salpeter IMF, where more low-mass stars than high-mass stars 
are formed, we have shown that the mass of the ISM gets ex- 
hausted more rapidly than for the IMFs favouring higher masses 
(see Sect. I4.31 l. This results in a steeper decline in the SFR, the 
SN rates and the SN dust injection rates. Thus, our models lead 
to lower dust masses than the models bv lDwek et al.l d2007l) for 
a Salpeter IMF. 

The model for a system with M ini = 1.3 x 10 12 M and M c \ 
= 800 M (Fig. 1 top left), a Larso n 1 IMF and the 'high ' SN 
efficiency is directly comparable to Valian teet alJ d2009l) . Our 
results are in good agreement with their results, although our 
model does not quite reach 10 8 M . This discrepancy may stem 
from differences in the SN dust yields u sed. Additiona l ly, the 
treatment of stars between 40-100 M bv lValiante et al.l (2009) 
is not unambiguously traceable. 



5. Discussion 

We have shown that, depending on the assumptions for certain 
model parameters, dust masses in excess of 10 8 M can be ob- 
tained in our model. All models leading to dust masses > 10 8 M 
within the first 400 Myr are listed in Table|4] In Fig.Q~3]we show 
the resulting relations between the dust mass and stellar mass, 
SFR, and metallicity at 400 Myr. 

It is evident that the dust yields are related to the mass of the 
galaxy. For a given combination of IMF, SN dust production ef- 
ficiency and M c i higher dust yields are obtained with increasing 
Mi n i. In the discussion below we consider thigh with M c \ = 100 
M as the reference case. 

The left panel in Fig. [13] depicts the dust mass versus the 
stellar mass. For a given IMF the stellar masses increase with 
Mi n i. However for a given M;„i, a wide spread of the stellar mass 
is found with varying IMF. We find that for the reference case, 
galaxies with Mj n ; > 5 x 10 11 M and an IMF biased towards 
higher masses reproduce dust masses > 10 8 M , while this is 
only possible in the most massive system with a Salpeter IMF. 
The stellar masses for these galaxies and with the more top- 
heavy IMFs are in good agreem ent with observations of galax- 
ies athigh-z. lSanfini et all ( 1201 Oh and lMichalowski et alJ d2010a) 
shows that, out of a sample of high-z submillimetre galaxies and 
local ULIRGs, most of these are assembled at stellar masses 
around 10 11 M and have dust masses of 10 8 ~ ll) M . 

Alternatively, a 'maximum' SN dust production efficiency 
can account for the required dust mass. For this SN efficiency, 
even a higher dust destruction in the ISM of M a = 800 M 
can be accomodated. However, assumptions of either e max or 
M c i = 800 M are both controversial. So far only in the SN 
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Fig. 11. Evolution of the stellar mass and mass of metals based on EIT08M. Results are shown for M; n i = 1.3 x 10 12 M (first 
column), 5 x 10" M (second column), 1 x 10" M (third column), and 5 x 10 10 M (fourth column). Top row: Stellar mass 
M,(t). Bottom row: Mass of metals Mz(t). The red lines mark a metal mass of 10 8 M . The solid, dotted, dashed, dashed-dotted, 
and dashed-dot-dotted lines represent the Salpeter, mass-heavy, top-heavy, Larson 1, and Larson 2 IMF, respectively. 



remnants Cas A (IWilson & Batrlall2005t iDunne et al.l l2009) and 



Kepler (Gome z et al. 20091) havedust masses been claimed that 



are consistent with 'maximum' SN efficiency, e max . Most SN 
dust observations reveal dust masses implying efficiencies in the 
range of ei nw -ehi?h (see G AH1 1). Additionally, theoretical mod- 
els ( Nozawa et al.l 12010) predict that depending on the density 
and geometry of the CSM and shocks, only a part of the dust 
may survive the destructive reverse shock. 

The uncertainty of dust destruction through SN shocks has 
been discussed in Sect. 13.51 The results obtained in this work 
strengthen the concerns of a high value for M c \. Apart from the 
'maximum' SN efficiency, for a high value of M c \ = 800 M , 
only the combination of a 'high' SN efficiency thigh with M m \ = 
1.3 x 10 12 M leads to dust masses close to 10 8 M , regardless of 
the IMF. However the remaining gas mass Mc(t) is > 10 12 M , 
implying dust-to-gas mass ratios ?7d(f) < 10~ 4 , is not supported 
by observations. The low dust-to-gass mass ratio is also, why 
dust destruction in the ISM is less efficient than in the lower 
mass galaxy systems. 

In principle, the fairly sensitive interplay between dust-to- 
gas ratio ?7d(f) and SN rates might be important when contem- 
plating additional dust source in less massive galaxies. As long 
as SN rates are high, any increase in the dust mass, and hence 
in rjd(t), leads to higher destruction rates. This shows that for 
M c i = 800 M it is difficult to reach high amounts of dust, un- 
less a rapid enrichment with high dust masses from either stellar 
sources or grain growth in the ISM takes place. In the latter case, 
grain growth rates must be comparable to SN dust injection rates 
for e max . Alternatively, a lowering of the dust-to-gas ratio rj^it) 
due to, for example, infalling gas might be an option. However 
this alternative results in an increased total mass of the galaxy 
anyway. 

Given these uncertainties, we find that the reference case 
constitutes the most likely scenario. In Fig. [13] (middle panel) 
we have plotted the dust mass versus SFR. The highest SFRs 
can be sustained for the high-mass weighted IMFs in the more 



massive galaxies. For the less massive systems, also with top- 
heavy IMFs, the SFR declines rather fast. In comparison to ob- 
servations of galaxies containing high dust masses at z > 6, we 
find that all systems more massive than 1 x 10 11 M can be 
broug ht into agreement with the SFRs derived from observation s 
(e.g. JBertoldieLani2003l:lDwek et al.ll2007tlWang et aD[2010h . 
The SFRs inferred from the observed far-infrared luminosities 
are naturally uncertain and sensitively dependent on assumptions 
about the IMF, the duration of the starburst, and the contributio n 
of the active galactic nucleus (AGN) (e.g.. lOmont et al.ll2.Q01l) . 
Commonly a Salpeter IMF is us ed, leading to high SFRs up to 
a few times 10 3 M yr" 1 (e.g.. IWalter et all 12004: Wan g et al.l 
120101) . As pointed out bv lDwek etalJ d2007l) . the inferred SFR 
however decreases to about 400 M yr~' if a top-heavy IMF is 
assumed. In Fig. [T3] (middle panel) it is evident that galaxies 
between 1-5 x 10 11 M G , and top-heavy IMFs result in SFRs be- 
tween about 300-700 M G yr" 1 . 

The right panel of Fig. Q~3] shows that the metallicity de- 
creases with increasing dust mass. IMFs favouring low masses 
lead to lower metallicities than the more top-heavy IMFs. The 
IMFs favouring low masses lock most material used for star for- 
mation in low-mass stars. These stars, however, do not recycle 
material back into the ISM, which therefore gets rapidly ex- 
hausted. In turn this leads to a high stellar mass M*(t). The more 
top-heavy IMFs form more short lived massive stars at higher 
rates. These stars recycle a copious amount of their mass back 
into the ISM. Therefore these IMFs lead to lower stellar masses 
M*(f), since the remnant mass of SNe is low. In addition, mas- 
sive stars enrich the ISM with metals, while owing to their longer 
lifetimes, low-mass stars do not eject a significant amount, if any 
at all, of heavy elements back into the ISM within 1 Gyr. This 
leads to lower metallicity for the lower mass star weighted IMFs, 
and higher metallicity for the top-heavy IMFs. 

Metallicities > Z have been found in strong star-forming 
galaxies, such as ULIR Gs or submillimetre galaxies, as well 
as in high-z QSOs (e.g.. iFan et aDl2003t iFreudUng et al.ll2003b 
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Fig. 12. Comparison of the evolution of the dust mass between 
EIT08M including and without the SMBH formation. The initial 
gas mass of the galaxy M m \ = 5 x 10 I() M . Calculations are 
shown for a 'maximum' (top row) and a 'high' (bottom row) SN 
dust production efficiency e\{m). Dust destruction is taken into 
account for M c \ =100 M (left column) and M Q (right column). 
Calculations including the SMBH are performed for a SMBH 
mass Msmbh = 5 x 10 9 M (solid curves) and for Msmbh = 3 x 
10 9 M (dotted curves). The calculations without the SMBH are 
illustrated as dot-dashed curves. The black lines represent the 
Salpeter IMF and the magenta lines the Larson 2 IMF. The grey 
horizontal dashed line marks the limit of 10 8 M of dust. The 
grey vertical dashed lines indicate epochs at 100, 400, and 500 
Myr after the onset of starburst. 



iKawara et alj |2010|) . In comparison to our calculated models, 
such metallicities can only be reached with the IMFs biased 
towards higher masses. In the lowest mass systems with these 
IMFs, the metallicity reaches values of ~ 5 x Z . 

5. 1. Caveats in our approach 

The presented models are based on rather simple assumptions, 
such as a closed box environment, a common constant initial 
SFR of the starburst, and a very simple treatment of the SMBH 
growth. With these assumptions the models are comparab le to 
similar works (e.g.. lDwek et alj[2007l : IValiante et alj |2009). and 
diverse evolutionary trends could be investigated in more detail, 
as discussed in previous sections. 

However, such models usually do not capture possible im- 
pacts on the evolution of the galaxy arising from galaxy merg- 
ers or gas flows powered by, for example, SNe or the SMBH. 
Different evolutionary paths likely caused by these effects and 
the mass of the galaxy possibly result in starburst intensities dif- 
ferent from our assumption, in turn leading to a different tempo- 
ral progression of various quantities. Furthermore, the growth of 
the SMBH is linked to several physical processes, e.g., the en- 
ergy feedback of SNe, and realistically may not take place at a 
constant growth rate. Despite the neglect of these effects in our 
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model, some general remarks about their influence can be made 
based on our results. 

For example, infall rates that could affect the systems might 
have to be unrealistically high because of the very high SFRs. 
However, merging galaxies might have an effect. The same ap- 
plies to outflows in the host galaxy or the feedback from the 
SMBH into the ISM. In our models the loss of material to fuel 
the SMBH can in principle be interpreted as some constant 'out- 
flow' of the host galaxy. Although we introduced a rough treat- 
ment for the SMBH growth, we could show that 'outflow' rates 
of the host galaxy of about 7.5-12.5 M yr~' on a timescale 
of 400 Myr only affect the amount of dust in the least massive 
systems (and IMFs biased towards low-mass stars), but are neg- 
ligible in the larger galaxies. In this regard, it remains to be in- 
vestigated whether the energy deposition by SNe/AGNs in the 
ISM could be high enough to initiate sufficiently higher outflow 
rates to impact the total amount of dust in a galaxy. 

A change in the evolution of the dust mass and other prop- 
erties may also result from quas ar winds, in whi ch dust for- 
m ation has also bee n suggested (Elvis et al. 1 12001 . According 
to lElvis et akl d2002l) . the estimated mass-loss rates of about > 



10 M yr 
10 47 ergs s~ ' 



in the most luminous quasars with luminosities > 
(e.g.. lOmont et al.ll200ll l200l ICarilli et alJl200U 



lBertoldi & Cox 2002) imply an amount of ~ 10 7 M of dust pro- 
duced over 10 8 yr. Comparing these rates to our model results, 
we find that neither the mass-loss rates nor the dust mass seem 
to be high enough to be relevant for the evolution of dust and 
other properties. 

We find that a strong impact on the evolution of dust is 
caused by the dust sources. Despite the included detailed treat- 
ment of the dust contribution from stellar sources in the devel- 
oped model, the poorly understood dust production by SNe, but 
also dust destruction by SN shock interactions, very likely con- 
stitute the largest uncertainties in the evolution of dust. Although 
not included in the model, alternative dust sources su ch as dust 
grain growth in the ISM might be relevant (e.g., iDwek et all 
l2007HDrainell2009HMichalowski et alJl2010bh . 
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Fig. 13. Relation between dust mass, stellar mass, SFR, and metallicity in correlation with the total mass of the galaxy and the 
IMF for an epoch at 400 Myr. Left panel: Dust mass M&(t) versus stellar mass M,(t). Middle panel: Dust mass Md(f) versus SFR 
if/(t). Right panel: Dust mass M^(t) versus metallicity Z(f). The dark grey shaded area marks the dust masses M^{t) obtained with 
a 'maximum' SN dust production efficiency e max (m) and no dust destruction in the ISM {M c \ = M ). The light grey shaded area 
marks dust masses M^{t) for the 'low' SN dust production efficiency ei ow (m) with dust destruction in the ISM (M c i = 800 M ). 
Cross points of the thin white lines indicate the data points for each IMF and mass M; n ; of the galaxy. The coloured filled circles 
signify the dust masses obtained for a 'high' SN efficiency ehigh('w) and a dust destruction in the ISM with M c \ = 100 M . The size 
of the circles is scaled by the initial mass M m \ of the galaxy. The black, green, cyan, magenta, and blue colors denote the Salpeter, 
mass-heavy, top-heavy, Larson 1, and Larson 2 IMF, respectively. 



6. Conclusions 

In this work we have developed a chemical evolution model for 
starburst galaxies at high redshift. The main purpose is to in- 
vestigate the evolution of the dust content arising from SNe and 
AGB stars on timescales less than 1 Gyr. In addition, we elabo- 
rated on the evolution of several physical properties of galaxies. 
The model allows the exploration of a wide range of parame- 
ters. The main parameters varied in this work are the mass of the 
galaxy, the IMF, SN dust production efficiencies, and dust de- 
struction in the ISM through SN shocks and stellar yields. The 
main results can be summarized as follows. 

1. The total amount of dust and the physical properties of a 
galaxy are strongly dependent on the IMF and correlate with 
the mass of the galaxy. For many properties we find an in- 
creasing disparity between the IMFs with decreasing mass 
of the galaxy. Higher dust masses are obtained with increas- 
ing mass of the galaxy. 

2. The maximum dust masses can be obtained with IMFs bi- 
ased towards higher mass stars and a 'maximum' SN dust 
production efficiency in galaxies with masses between 5 x 
10 10 M and 1.3 x 10 12 M Q . This is found independently of 
the strength of dust destruction in the ISM. The case of 'max- 
imum' SN efficiency and no destruction constitute the max- 
imum possible dust masses attainable with SNe and AGB 
dust production. The case with 'low' SN efficiency and a 
high dust destruction (M c i = 800 M ) gives the lowest possi- 
ble dust masses. 

3. The contribution by AGB stars is most visible in cases where 
a 'low' SN dust formation efficiency is considered. In this 
case the mass-heavy IMF dominates in all considered initial 
masses M; n ;, and leads to the highest dust masses. However, 
dust production with this efficiency is found to be insuffi- 
cient to fully account for dust masses in excess of 10 8 M 
within 400 Myr. In the galaxies with M ini > 5 x 10 11 M Q , 
this limit can be reached at epochs > 400 Myr if there is no 
dust destruction in the ISM. 



4. The total dust mass is considerably reduced when destruction 
by SN shock waves is taken into account. The strength of 
destruction in the ISM is given in terms of the mass M c \, 
which is the mass of the ISM that is swept up and cleared 
of the containing dust through one SN remnant. Top-heavy 
IMFs are sensitive to dust destruction, while IMFs favouring 
lower masses are more resistent to destruction due to lower 
SN rates. This leads to a significantly greater reduction and, 
in some cases, to a lower amount of dust with the top-heavy 
IMFs. 

5. At early epochs (< 200 Myr) SNe are primarily responsible 
for a significant enrichment by dust. For a 'high' SN dust 
production efficiency SNe can generate 10 8 M within the 
first 100 Myr. In comparison, for 'low' SN dust production 
efficiency, SNe are only able to increase the dust mass up to 
a few times 10 6 M in the first 100-150 Myr. 

6. Taking the growth of the SMBH into account leads to a re- 
duction of the amount of dust of at most ~ 50 %. This is 
achieved for a SMBH mass Msmbh = 5 x 10 9 M in a galaxy 
with M mi =5 x 10 10 M and only for IMFs favouring low- 
mass stars. In more massive systems and for top-heavy IMFs 
no variation in the dust mass and properties of the galaxy is 
encountered. 

7. To account for dust masses > 10 8 M we find that galaxies 
with an initial gas mass of M m i = 1-5 x 10 11 M in con- 
nection with top-heavy IMFs, 'high' SN dust production ef- 
ficiency ehigh, and dust destruction in the ISM of M c i = 100 
M are favoured. Models with the highest amount of M m i = 
1.3 x 10 12 M and lowest of M ini = 5 x 10 10 M are discour- 
aged. 

More refined models including diverse flow scenarios tracing 
self-consistent SF histories, detailed SMBH growth scenarios, 
and several plausible dust sources are major future develop- 
ments. 
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